Integrability of the double pendulum — the Ramis' question 
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In this short note we address the problem of integrability of a double pendulum in the con- 
stant gravity field. We show its non-integrability using the combination of algebraic and numerical 
approaches, namely we compute the non-commuting generators of the monodromy group along a 
particular solution obtained numerically. 
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I. INTRODUCTION 



A double pendulum is a system that is often modeled as two mass points connected between themselves by weightless 
inextensible rods. Mathematically this corresponds to considering the motion of two points ri and r 2 subject to the 
following constraints: 



pa = (r x - r ) 2 - l\ = 0, ^ 2 = (r 2 - ri ) 2 - l\ = 0, 



the point Tq being fixed. 




FIG. 1: Double pendulum - angle parametrization 



A convenient approach for the integrability analysis would be resolving this constraints, i.e. introducing the gener- 
alized coordinates automatically taking the constraints into account. As it can be seen from the figure [1] the system 
can be parametrized by two angles. That is for Liouville integrability one needs two independent first integrals. The 
free flat motion of a double pendulum is a classical example of an integrable system: the corresponding first integrals 
are the energy and the component of the angular momentum orthogonal to the plane of the motion. But already in 
the presence of the gravity field the situation is less clear. It has been shown (mainly using the techniques like the 
Poincare sections) that the system might possess some chaotic behaviour. The question of Jean-Pierre Ramis sounds 
roughly as: "Can one rigorously prove non-integrability?" 

In this note we discuss meromorphic non-integrability of the double pendulum within the framework of the algebraic 
approach inspired by the results related to analysis of the monodromy group (Q) and the differential Galois group 
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). As a mechanical system the double pendulum in the constant gravity field g can be described by the following 
Lagrange function: 

L = a\ + a\a.2 cos(«i — a 2 ) + + 2<?cos(ai) + .gcos(a 2 )- (1) 

Here and in what follows we assume all the masses as well as all the lengths of the constraints to be equal to 1, this is 
done to simplify the formula, but certainly the analysis can be carried out with arbitrary values of these parameters. 



II. DESCRIPTION OF THE METHOD 



The standard procedure of application of the methods of Ziglin or Morales-Ramis for recovering the obstructions 
to integrability consists of complexifying the equations and analyzing the monodromy group or respectively, the 
differential Galois group, of the system of variational equations along an explicitly known particular solution of the 
initial dynamical system. The idea is that for an integrable system the obtained group should not be too complicated, 
i.e. up to some very special cases should be (virtually) abelian (see [l[ for details and examples of application). 

The key difficulty is that the particular solution should on one hand be rather "nice" so that one would be able 
to compute the corresponding monodromy or differential Galois group, and on the other hand not trivial so that the 
computed group would contain a sufficient amount of sources of non-commutativity. For example for the system of 
Euler-Lagrange equations governed by the Lagrangian L (eq. [1]), that interests us, there are no obvious particular 
solutions except the equilibria. 

To handle this difficulty we have suggested (0) a method of computing the generators of the monodromy group 
based on a particular solution obtained numerically. It consists of basically searching for the paths in the complex 
plane that correspond to finite order branching points of the solution of the initial system. Having found several of 
them we solve the system of variational equations in parallel with reconstructing the Riemann surface of the particular 
solution and obtain the corresponding monodromy matrices. If among them we see non-commuting ones it permits 
to conclude meromorphic non-intcgrability. 



III. RESULTS OF THE COMPUTATION 



We are thus solving the system of differential equations 

x = v(x), S = A(x)H, (2) 

where " ' " denotes the derivative with respect to the complex time. The first equation defines the evolution of 
x = (ai, Q!2, Ax, «2) T with the right-hand-sides obtained from the Euler-Lagrange equations governed by L (equation 
[!}; the second one is the matrix system of variational equations with A(x) = 

Starting from {ot\ , q?2, Qt\ , 0*2) = (0.1,-0.3,0.2,0.4) and following the paths 71 and 72 encircling the points t\ = 
0.5 + 0.9z and ti — 0.5 — 0.9* (three times each of them - see fig. [5]) one obtains respectively the monodromy matrices 
Mi and M 2 . 

14.55+ 10.90* \ 
12.91 - 7.99i 
-8.18-5.42i 
8.59 + 9. 31* / 



/ 20.72- 17.12* 
-17.67+ 12.78* 
12.28 + 6. 91* 
\ -11.84- 12.56* 



15.79- 1.34* 
-12.24-0.06* 

3.55 + 7.78i 
-1.31 - 10.39* 



4.94 + 29.46? 
-1.93-24.87* 
-13.07 + 7.86? 

19.32 - 4.06* 




-18.72- 17.12* -15.79- 1.34* -4.94413 + 29.46i 14.55 + 10.90* \ 

17.67+ 12. 78i 14.24 - 0.06i 1.92944-24.87* -12.91-7.99* 1 

-12.28 + 6.91* -3.55 + 7.78i 15.0698 + 7.86*' 8.18-5.42* 

11.84- 12.56* 1.31 -10.39* -19.319-4.06* -6.59 + 9.31*/ 



One easily checks that these matrices do not commute, that shows meromorphic integrability. 
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FIG. 2: Loops in the complex plane corresponding to the obtained monodromy matrices 

It is worth noting that all the computations are done with a controlled precision, so the result is certainly more than 
just a numerical evidence. One should also make a remark that the computation passes some obvious compatibility 
tests. For instance if one continues integrating ([2]) along 71 or 72 for three more loops, one obtains the corresponding 
matrices squared. The eigenvalues are also coherent with their symplectic nature. Moreover one can notice, that the 
obtained matrices can be written in the form 

Mi = Id + A + iB, M 2 = Id- A + iB, 

for some real matrices A and B, and their commutator [Mi,M2] — 2i[A,B]. This observation also simplifies the 
precision control. 



Conclusion 



As we have seen in this note the (computer assisted) proof of non-integrability of the system describing the motion 
of a double pendulum is possible. It would be interesting to see if the difficulty of the procedure depends on parameters 
of the system and in particular if it is related to the value of the energy first integral. 
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